Towards a formalism for mapping the spacetimes of massive compact objects: 

Bumpy black holes and their orbits 

Nathan A. Collins and Scott A. Hughes 
Department of Physics, MIT, 77 Massachusetts Ave., Cambridge, MA 02139 

Astronomical observations have established that extremely compact, massive objects are common 
in the universe. It is generally accepted that these objects are, in all likelihood, black holes. As 
observational technology has improved, it has become possible to test this hypothesis in ever greater 
detail. In particular, it is or will be possible to measure the properties of orbits deep in the strong field 
of a black hole candidate (using x-ray timing or future gravitational- wave measurements) and to test 
whether they have the characteristics of black hole orbits in general relativity. Past work has shown 
that, in principle, such measurements can be used to map the spacetime of a massive compact object, 
testing in particular whether the object's multipolar structure satisfies the rather strict constraints 
imposed by the black hole hypothesis. Performing such a test in practice requires that we be able 
to compare against objects with the "wrong" multipole structure. In this paper, we present tools 
for constructing the spacetimes of humpy black holes: objects that are almost black holes, but that 
have some multipoles with the wrong value. In this first analysis, we focus on objects with no 
angular momentum. Generalization to bumpy Kerr black holes should be straightforward, albeit 
labor intensive. Our construction has two particularly desirable properties. First, the spacetimes 
which we present are good deep into the strong field of the object — we do not use a "large r" 
expansion (except to make contact with weak field intuition). Second, our spacetimes reduce to the 
exact black hole spacetimes of general relativity in a natural way, by dialing the "bumpiness" of the 
black hole to zero. We propose that bumpy black holes can be used as the foundation for a null 
experiment: if black hole candidates are indeed the black holes of general relativity, their bumpiness 
should be zero. By comparing the properties of orbits in a bumpy spacetime with those measured 
by an astrophysical source, observations should be able to test this hypothesis, stringently testing 
whether they are in fact the black holes of general relativity. 

PACS numbers: 04.25.Nx, 04.30.Db, 04.70.Bw 

I. INTRODUCTION 
A. Motivation 

Observations have now established that the cores of nearly all nearby galaxies contain a massive, compact, dark 
object These objects range in mass from several IQ^ Mq to several 10® A/q. Extremely compact stellar mass 

objects (~ \QMq or so) exist and have been studied in the galactic field (see, e.g., Ref. |3] for a review). Evidence 
suggests the existence of objects with intermediate masses, 10^ Mq — 10'' Mq, filling the gap between the supermassive 
and stellar mass objects (see, e.g., Ref. for a review). The most generally accepted explanation is that these compact 
bodies are massive black holes. 

Although this is the most generally accepted explanation for these objects, it is not the only explanation. In 
some cases, the massive dark objects seen in galaxy cores can be explained quite well as dense clusters of stars or 
stellar remnants. Such models are rapidly becoming disfavored in many cases as our ability to study the central 
regions of galaxies improves — many of these putative clusters would have to be so compact that they would not 
be gravitationally stable. By using "exotic" matter, it becomes possible to build objects that are massive, compact, 
but stable. For examine, by tuning the mass and self interaction of a massive scalar field one can build an 

object that is consistent with much of the observational evidence available today. Indeed, the fields that describe some 
of these black hole alternatives are similar to some dark matter candidates, leading to the suggestion that massive 
compact objects could be dark matter condensates rather than black holes. 

Other recently proposed black hole alternatives are motivated by a desire to avoid the information paradox — the 
loss of information through the black hole's event horizon. Such models find ways of eliminating the event horizon 
altogether, for instance by replacing the event horizon with a hard surface surrounding a ball of negative energy density 
(the "gravastar" model) 0> or by postulating that spacetime itself undergoes a phase transition in the presence of 
very strong gravitational fields 0|. If such objects exist in nature, they should have a deep, strong field structure very 
different from that of black holes |^ . 

Astronomical measurements are now becoming able to probe into the very strong field of compact objects: optical 
and infrared observations track stellar orbits at the core of the Milky Way, probing the spacetime of the presumed black 
hole at Sgr A* ,11, .12J ; x-ray observations of quasi-periodic oscillations from black hole candidates carry information 
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about gas in the hole's deep strong field [13|; and future gravitational- wave observations ma y b e able to track the 
sequence of orbits followed by a compact body that slowly spirals into a massive black hole 0, 0| . The question 
of whether these objects are truly the black holes of general relativity or are described by some alternative model 
reduces to the question of how one may use measurements of orbital properties to map the spacetime structure (i.e., 
the gravitational field) of these objects One thus needs to be able to relate the properties of the measured orbits 
to the structure of the central gravitating objects. A powerful way of doing this is by a multipole expansion of the 
compact object's spacetime. 



B. Multipoles of massive compact objects 

In Newtonian theory, the gravitational field of a body is simply described by expanding the potential in spherical 
harmonics. The potential $ must satisfy 

V^'l' = inGp (interior) , 

~ (exterior) . (1-1) 

In the exterior, the potential may be expanded as 

* = -GE^^- (1-2) 

Im 

By matching to an expansion of the interior solution and enforcing Eq. , we see that the coefficients M;™ are mass 
multipole moments: numbers that describe the angular distribution of matter inside the star. For simplicity, let us 
focus for a moment on axially symmetric objects, so that only m — matters. Then, for example, Mqq = M, the total 
mass of the object. By an appropriate choice of the center of our coordinate system, we put the moment Miq = 0. The 
first interesting moment is M20, the quadrupole moment of the object. This moment has the form QML^, where L 
is the object's "size" (e.g., its mean radius), and the dimensionless number Q describes the quadrupolar deformation. 
Higher moments can likewise be interpreted as ^-polar moments of the mass distribution. Because these moments 
directly determine the gravitational potential outside of the gravitating object, one can measure properties of its mass 
distribution by measuring the "shape" of the gravitational potential. Doing so by studying the properties of satellite 
orbits is the science of geodesy. 

Somewhat remarkably, such a description describes the exterior spacetimes of bodies in general relativity as well. For 
any gravitating body that is stationary, axisymmetric, and reflection symmetric across the equator — encompassing 
black holes plus a wide variety of perturbations — the exterior spacetime is fully specified by a pair of multipole 
moment families: the mass multipole moment Mi, plus the current multipole moment Si |l7l Il8l|. (Since we have 
restricted ourselves to axisymmetry, we only consider m = 0. We henceforth drop this subscript.) The current 
moment Si does not appear in Newtonian theory; it refiects the fact that mass and energy flows gravitate in general 
relativity. For example, the moment Si is the magnitude of the spin angular momentum of the body. 

If the gravitating body is a Kerr black hole, then the values of the mass and current moments are strongly restricted: 
in units in which G = 1, c = 1 (which we use throughout this paper), we must have (18j 

Mi+iSi^ M{iay . (1.3) 

This equation tells us that Mq — M, the total mass of the Kerr black hole, and Si — aM, the magnitude of the 
black hole's spin angular momentum in these units — precisely what we already expect. More interestingly, all higher 
moments are completely determined by these two values. The exterior spacetime of a Kerr black hole is completely 
determined by its two lowest multipole moments — its mass and spin. 

This is nothing more than a restatement of the "no hair" theorem |23, |^ H^l ■ By analogy with geodesy, this 
suggests that one can test the no hair theorem by measuring orbits near black holes. Using a spacetime that does 
not necessarily assume the Kerr form of the moments, one could then determine Mi and Si. If that object is in fact 
a black hole as described by general relativity, the only free moments are those for / — and I — 1. Once they have 
been determined, all higher moments are constrained via Eq. 1)1. Such "geodesy for black holes" (which has been 
given the names "bothrodesy" and "holiodesy" ^^) would provide a stringent test of the black hole nature of massive 
compact objects in the universe. 

The first detailed analysis of how one might be able to falsify the black hole nature of a massive compact object 
was by Fintan Ryan |24j . Ryan showed how to build the spacetime of an object with arbitrary multipole structure, 
and then analyzed the orbits of small bodies in that spacetime. ("Small" means that these bodies do not themselves 
significantly distort the spacetime, and so can be treated as following approximately geodesic trajectories). His analysis 
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demonstrated that the accumulated orbital phase was sensitive to these multipoles. Orbital phase (or some surrogate 
of this phase) is directly observed by, for example, x-ray timing (today) and gravitational- wave detectors (future). 
One could thus imagine using measurements of accumulated orbital phase to test the black hole nature of a massive 
compact object. Focusing on gravitational-wave measurements with the planned space-based laser interferometer 
LISA j25i, Ryan showed that enough multipoles should be measurable to easily falsify the black hole hypothesis. In 
many cases, enough multipoles would be measurable (up to ^ ^ 5 or 6) to stringently constrain the object's black hole 
nature p^ . 

Unfortunately, the multipole expansion used by Ryan does not work very well in the deep strong fields of massive 
black holes, where one expects orbital phases to most stringently test the black hole hypothesis. Multipole moments 
essentially label different powers of a 1/r expansion. In the strong field of a black hole (small r), such an expansion is 
not going to be very useful |27j . The in-utility of this expansion is reflected by the extremely large number of terms 
that must be kept to describe a spacetime with arbitrary multipole moments at small radius (cf. Ref. 0|, Sec. III). 

C. Bumpy black holes 

We advocate a different approach. The reason for introducing a multipolar expansion is to describe a candidate 
spacetime differing from that of a black hole. If one accepts as a starting point the idea that the black hole hypothesis 
probably describes the massive compact objects in question, then one just needs a spacetime to compare against that 
differs slightly from that of a black hole. Our goal is then to set up a null experiment: we find a trial spacetime that 
exhibits slight deviations from the spacetime of a black hole. If the black holes of nature are the black holes of general 
relativity, we will measure the deviation to be zero. 

Past work on candidate objects to test the black hole hypothesis has focused primarily on boson and soliton stars 
0,123. Though of great intrinsic interest, there is no natural way for a boson star spacetime to smoothly limit to 
the spacetime of a black hole. If the massive compact objects we observe in the universe are in fact black holes, then 
tests based on the boson star model will not provide useful constraints on orbit observations. As a "straw man" for 
the black hole hypothesis, boson stars may unfortunately contain too much straw. 

We advocate instead the use of bumpy black holes: objects that have a multipolar structure that is very nearly, 
but not quite, that of a black hole. As the name suggests, these are black holes with small bumps on them. If the 
universe's observed massive compact objects are in fact black holes, then we will find that the amplitude of these 
bumps is zero (within measurement uncertainty). A bumpy black hole should be a trial spacetime that behaves well 
deep into the strong field, and that exhibits a controllable deviation from the Kerr solution. In particular, these trial 
spacetimes should reduce to normal black holes when the the deviation is set to zero — bumpy black holes become 
normal black hole when the bumps are removed. This reduction to normal black holes is a crucial element of using 
bumpy black holes as a basis for a null experiment. 

A key piece of our guiding philosophy is that the notion of multipoles is most useful in the weak field of an object. 
One shouldn't be too attached to multipole moments if the goal is an analysis that applies to strong gravitational 
fields. Of course, by taking the weak field (large r) limit of the spacetimes we construct, bumpy black holes are very 
usefully described using multipoles. Indeed, the goal of our detailed calculations will be to assemble a perturbation 
that is purely quadrupolar when examined in the weak field. Our construction, however, works very well deep in the 
strong field, which is crucial for applying these notions to observations. 

An important question to ask at this point is how one can build a stationary spacetime corresponding to a bumpy 
black hole. A key portion of the proof of the no hair theorem demonstrates that any deformation to a black hole 
will tend to radiate very quickly, removing the bump and pushing us back to the Kerr black hole solution ,2L .23 | . 
Some mechanism must exist to maintain the bump. This is likely to require unphysical matter; the example which 
we describe in fact requires naked singularities. One might object that a bumpy black hole spacetime is thus, by 
construction, unphysical. Our viewpoint is that the physicality of these spacetimes is irrelevant. Our goal in this 
analysis is not to build a spacetime which might conceivably exist in nature. Instead, we wish to build a black hole 
straw man with just the right amount of straw to probe the nature of massive compact objects. 

D. Overview of this paper 

The goal of this paper is to present the bumpy black hole concept, to show how bumpy black hole spacetimes are 
generated, and to demonstrate that the magnitude of the bumps is encoded in the accumulated phase of the hole's 
orbits. We focus upon axisymmetric distortions of black holes — even in axisymmetry, an incorrect moment is enough 
to falsify the black hole hypothesis for a massive compact object. 
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We have argued that the language of multipoles is not useful for describing an object's strong field orbits. To 
substantiate this argument, we review the multipole description of spacetimes and their orbits in Sec.^ summarizing 
the key results of Ryan 24]. Ryan's formulae and the detailed description of the spacetime in the multipole language 
show that, as we try to characterize the massive object's strong field, the description becomes extremely complicated. 
Although it is possible in principle to use this description to develop tools for mapping spacetimes, it does not appear 
to be the best approach in practice. 

We then begin our detailed presentation of bumpy black holes. In Sec. IIIII we show how to build a bumpy black 
hole spacetime. The spacetime of a stationary, axisymmetric object is fully described by the Weyl metric 28], 

ds^ = -e^^'dt^ + e^^'-^-^idp" + dz^) + e-^^'p^dc^^ . (1.4) 

Our strategy is to pick an exact solution ip = ipQlfjz), 7 = 7o(/9, z) for which the line element (|1.4|) describes a 
black hole. For this first analysis, we specialize our background to Schwarzschild black holes; generalization to Kerr 
black holes should be straightforward in principle (though it may be somewhat involved algebraically). We then use 
this exact solution as a background against which to introduce a perturbation, putting ijj = + ipi and requiring 
ipi/i^o *C 1; a similar perturbation is defined for 7. The perturbations are constrained by the requirement that they 
satisfy the vacuum Einstein equations, expanded to first order. 

This formulation of the metric is particularly useful because the function tp reduces to the Newtonian gravitational 
potential of the source in the weak field. We therefore choose our perturbation tpi in such a manner that the weak-field 
perturbation can be thought of as changing the source's multipoles as measured in the weak field. We then solve the 
linearized Einstein equations in order to specify the perturbation throughout the exterior spacetime of the bumpy 
black hole. 

Before specifying our perturbations, we first examine the properties of orbits in the bumpy hole's equatorial (re- 
flection symmetry) plane fSec. Hvjl . Many useful quantities can be computed in terms of the perturbations V'l E^nd 
71 — the orbit's energy E, angular momentum L, and the location of the last stable orbit (a separatrix in orbital 
phase space dividing dynamically stable and unstable orbits). We also write down an equation describing the differ- 
ential advance of the orbit's periapsis. The periapsis shift arises from a mismatch between the radial and azimuthal 
orbital frequencies; as such, it can be a sensitive probe of the spacetime. Deviations of this shift from the canonical 
Schwarzschild value encode the black hole's bumpiness. 

We choose particular perturbations in Sees. IVland I VII A very simple and useful one is that of a point mass near 
the black hole. We build a bumpy black hole in Sec. [3 by placing a pair of point masses with mass m / 2 each near 
the hole's "north" and "south" poles. The same system was used by Suen, Price, and Redmount (SPR) l23] to set up 
an analysis of a black hole with a deformed event horizon. Our analysis is similar to that of SPR, though we do not 
focus on the region of spacetime near the horizon. We build a second type of bumpy black hole in Sec. IVII bv placing 
a ring of mass fj, about the hole's equator. 

As the analysis of Sees. IVl and I VII shows, both the polar point mass and the equatorial ring do indeed change the 
metric's quadrupole moment. We demonstrate this by calculating weak-field periapsis precession in these spacetimes 
and showing that the shift contains a term which is exactly what we expect for weak-field quadrupolar distortions 
(computed in Appendix . Unfortunately, these perturbations also change the metric's monopole moment (i.e., its 
mass). Fortunately, we can build a purely quadrupolar distortion by combining negative mass polar points with a 
positive mass equatorial ring, or vice versa. (Bearing in mind that our goal is to build trial spacetimes for testing the 
black hole hypothesis, the unphysicality of a negative perturbing mass is not a concern.) In Sec. IVIl"^ we show that 
the weak-field periapsis precession with this combined mass distribution is identical to that of a Schwarzschild black 
hole plus a quadrupolar deformation. The points + ring perturbation to a Schwarzschild black hole thus perfectly 
matches our requirements for a bumpy black hole. We investigate the strong-field character of this spacetime in Sec. 
IVllBI showing in particular that the hole's bumpiness is usefully encoded in the strong-field periapsis precession. 

Concluding discussion is given in Sec. IVIIII In particular, we outline further work that should be done to connect 
the bumpy black hole concept to future astrophysical observations. Chief among the tasks needed is a generalization 
to Kerr black hole backgrounds; some steps in this direction are outlined in Appendix IbI 

II. METRICS AND MULTIPOLE MOMENTS: AN OVERVIEW 

As stated in the Introduction, one can build a spacetime by specifying a set of mass and current multipole moments 
(Ml, Si). In actuality, one builds a spacetime from a set of coefficients ajk which determine orbital characteristics; 
the multipole moments can then be extracted from these coefficients. In this section we briefiy describe the details of 
this analysis, and discuss why we believe this is not the most effective way to map black hole spacetimes in practice. 
Much of our presentation is essentially a synopsis of Fintan Ryan's calculation; see Ref. |2^ for detailed discussion. 
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To begin, we must pick a spacetime sufSciently general to encompass the stationary, axisymmetric sources we wish 
to describe. Ryan begins with the Hne element 



ds^ = -F {dt - u d4>f + - [e^T (V + dz"^) + p2d02j 

F 



(2.1) 



The functions 7, and w depend on p and \z\. The radial coordinate p labels displacement from the source's 
symmetry axis; z labels displacement above or below the source's "equatorial" plane. By construction, the spacetime 
is stationary and axisymmetric {F , 7, and lo do not depend on t or 4>), and is reflection symmetric about the equatorial 
plane (dependence on |z|). 

For an axial and reflection symmetric spacetime, the metric functions F, 7, and lu can be generated from the Ernst 
potential £ . This function and a related complex function ^ are deflned via 



where the function is related to w by 



£ ^F + i^ 



uj{p,z) 



+ z2 - g 

p' 5* 



F2 dz 



dp' 



(2.2) 



(2.3) 



The quantity in parentheses under the integral is evaluated at constant z. The function 7 can be determined once F 
and uj are known: the vacuum Einstein equations tell us 



97 
J~p 

dz 



4F2 



1 p dFdF 
2F^~dp~dI 



dz 

1 F2 duj duj 

2 p dp dz ' 



4 p 



duj 
'dp 



duj 
dz 



(2.4) 
(2.5) 



Then, 7 can be found at any point in spacetime with an appropriate line integral. See Ref. j3y|, pp. 165 - 167 for 
further discussion. For our purposes, the main thing to note is that knowledge oi £ — or equivalently ^ — specifies 
the entire spacetime metric. 

The function ^ can be expanded as 



p^z^ 



2\0 + k 



(2.6) 



The index j is strictly even. If k is even, ajk is real; if k is odd, aju is imaginary. From these coefficients, it is 
relatively straightforward to extract the spacetime's multipole moments {Mi, Si) using an algorithm developed by 
Fodor, Hoenselaers, and Perjes (FHP) [sj- By a recursive procedure that involves repeatedly differentiating £ and 
gathering terms, FHP show that ajk can be written in terms of the mass and current moments [see Ref. |2J|, Eqs. 
(35) - (41)]. For example, one can show that 



ao; = Ml + iSi + LOM , 

aio = (-1)'/2^^^M,+L0M, 

az_i,i = i(-l)('-i)/2_i!L_^,+L0M 



(2.7) 
(2.8) 

(2.9) 



"LOM" is an abbreviation for "Lower Order Moments" : a complicated (but known) sum of various combinations of 



Mj and S, 



Ryan [2J| uses the coefficients ajk and their relationship to (A/;, Si) as the basis for his spacetime mapping procedure. 
Since the spacetime is fully determined by ajk, it follows that its orbits are likewise determined. An orbit in the 
spacetime (|2.1f) is governed by three orbital frequencies: Sl^, related to the time to cover 27r radians of azimuth; ilp, 
characterizing oscillations in the p coordinate; and f2z, for oscillations in the z coordinate. Ryan shows (via a power 
law expansion using orbital speed as an expansion parameter) how these frequencies depend on ajk- Measurement 
of (ri^, rip, Jlj.) (or accumulated phases associated with these frequencies) can therefore be used to measure ajk- 
Ryan then inserts the measured values of ajk into the FHP algorithm, determining the multipole moments of the 



with j < I and k < 
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spacetime. Using this procedure, Ryan showed that future gravitational-wave measurements with LISA should be 
able to determine enough moments to strongly constrain the black hole nature of massive compact objects j26l] . 

Though sufficient to prove the principle, we believe that this procedure is not useful in practice for mapping 
the spacetimes of objects believed to be black holes. Referring to Eq. (|2.t)|l . we see that the coefficients ajk are 
essentially labels for an expansion in inverse distance. For strong field orbits (p, z ~ a few x M), a large number of 
these coefficients must be kept in order to model the spacetime with sufficient accuracy. One might hope that the 
coefficients ajk become small for large j and fc, making it possible to truncate the expansion of ^. This is not the case: 
because of the coupling to lower order moments [cf. Eqs. (|2.7|l - H2.9|l ]. these coefficients generically remain large even 
if the body has only a small number of non-zero multipole moments. 

In this language, the description of the spacetime and hence of its orbits becomes extremely complicated in the 
strong field. This makes testing whether a spacetime is close to that of a black hole very difficult. Naively, one might 
imagine requiring that a spacetime have the multipole moment structure 

Ml + iSi = M{iay + 6Mi + iSSi , (2.10) 

and then developing a formalism similar to that described here to place observational limits on the deviations 
{SMi,5Si). (Indeed, this how we originally conceived of this analysis.) One quickly discovers that the algebraic 
complexity associated with this approach is immense. Though no doubt possible in principle, a multipolar prescrip- 
tion like Eq. H2.10|) does not easily translate into a practical scheme for constraining the properties of massive compact 
objects. 

The lesson appears to be that multipoles, though conceptually clean and offering a beautiful description of weak-field 
gravity, simply aren't the best tools in the strong field. This shouldn't be surprising in a nonlinear theory like general 
relativity: since multipoles are basically labels in an inverse distance expansion (as we have repeatedly emphasized), 
a description that is clean in the weak field can easily turn into a mess when the nonlinearities are important. 



III. BUILDING THE SPACETIME OF A BUMPY BLACK HOLE 



Keeping in mind that the goal is just to build some candidate spacetime to be used as a straw man in testing the 
black hole hypothesis, we advocate a different approach. Our goal is to develop a family of spacetimes corresponding 
to stationary perturbed black holes: bumpy black holes. These spacetimes include black holes as a subset — we 
simply set the magnitude of the perturbation to zero. We construct these spacetimes in a manner that makes an 
exploration of its strong-field properties simple. 

We will focus on axisymmetric spacetimes, in keeping with our goal of analyzing axisymmetric deformations of 
Kerr black holes. For this first analysis, we will further specialize to deformations of Schwarzschild black holes. 
Stationary, axisymmetric deformations of Schwarzschild black holes were in fact studied by SPR 29] with the goal 
of characterizing distortions to the event horizon. (Their analysis was a part of the "Membrane Paradigm" ; sec Ref. 
[321.) Their calculation makes an ideal starting point for our analysis; the following discussion closely follows Ref. 
Ha, Sec. IIIA. _ 

As mentioned in the Introduction, the spacetimes we consider are described by the Weyl metric (28l |: 



(3.1) 



The vacuum Einstein equations for this line element reduce to 



dp 

d-f 
dz 



+ 



1 dij 
P dp 



9-0 



dip dip 
dp dz 



dz^ 

2 



0, 

&ip_ 
dz 



(3.2) 
(3.3) 
(3.4) 



Equations and (|?^ are identical to Eqs. and with w = and = e^'^. Note that E^^ implies 

^ is a harmonic function in a fictitious Euclidean space with cylindrical coordinates p, z, and (p |29ll3Cl| . 

Following SPR we observe that the function ip is the Newtonian potential far from a source. A reasonable way 
to perturb tp is by adding potential terms that correspond to particular mass distributions perturbing the background 
black hole; distributions that change the weak-field multipole structure of the hole are particularly interesting. [Note 
the linearity of Eq. H3.2|) — exact solutions for ip can be constructed by superposition.] We only require that the mass 
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/i of the perturbations be small compared to the mass M of the black hole, allowing us to expand and 7 to first 
order. 

Notice that, with our form of the perturbation, we expand in powers of /x/Af, rather than l/r (where r is some 
measure of distance from the source). The approximation we introduce should thus be well-behaved for any r, 
including into the strong field. Notice also that, because Eqs. (|3.2ll - (|3.4|l come from the vacuum Einstein equations, 
our metric will only hold where there is no matter. To perturb the background black hole, we will add matter in the 
form of point particles and one-dimensional rings. The metric will therefore not apply at the points containing that 
matter; it will in fact be singular at those locations. As long as we only examine regions of spacetime external to 
these perturbing sources, this singular behavior does not pose any difficulty. 

For describing black holes, it is convenient to use prolate spheroidal coordinates u and v. 

p — Msinhusint! , (3.5) 
z = Mcoshwcos?;. (3-6) 

The coordinate v € [0, tt] is a polar angle; u S [0, 00) is effectively a radial coordinate. These coordinates cover the 
entire exterior Schwarzschild spacetime: the coordinate u = maps to the event horizon, r — 2M [cf. Eq. (|3.15|) 
below]. In the Weyl coordinates (p, z), this corresponds to a cylindrical rod at p = running from z — ~M to z = AI. 
The line element l|3.1|l becomes 

ds^ = -e^'^dt^ -t-M2e2''"2^(sinh^u + sin2w)(dM2 +M2e-2'^sinh2usin2ud02 _ (^g j) 

We now put xjj ~ ipo + tpi, 7 = 70 +7i- We take the perturbations (V'i,7i) to be small compared to (ipo^jo)- 

V'l ~ ^V'o , (3.8) 
and likewise for 7. The Einstein equations constraining 7, H3.3|) and H3.4|l . become 

^71 _ ^ 2 „^ f dipo dipi , dipodipi\ ^ 9-00 9-01 d-ipo d-ipi 



— = 2p cotv [ — — + — — + 2pztant> -^-^ ^ ^ , (3.9) 

ou \ op oz oz op J \ Op Op Oz oz ' 

dii „ 2 _Y90o 901 i90o90i\ ^ „ /^9-0n9-0i , 90o 90i 



„ -2p'coiv[^^-^^\+2pzis.nv[^^ + ^^\. (3.10) 
ov \ Op Op oz oz J \ op Oz Oz op J 

We have only kept terms to leading order in the perturbation. It turns out that these two equations will overdetermine 
the solution (cf. |30j |. p. 167 for discussion). For our purposes, this means that we need integrate only one of them to 
determine 71 (up to a constant of integration). We will use Eq. (|3.10() : the solution we find also satisfies H3.9|l . as is 
easily verified by direct substitution. 

At this point, we specialize our background to the Schwarzschild metric: we put [2^ 

00 = lntanh-u/2 , (3.11) 
1 , / sm V 



70 = - In 1 + —^ . (3.12) 
^ V smli u / 

Using Eqs. H3.11|l and ()3.12f) and changing all instances of p, z to u,v as appropriate, Eq. (|3.1U|) becomes 

971 2 [tan V (501 /dv) — tanh u (90i / du)] 
dv sinhw (cothw tan u -|- tanh w cot v) 

For completeness, the equation for d"fi/du becomes 

971 2 [cot V (901 /du) -I- tanh u (90i /dv)] 



(3.13) 



dv sinhu (cothu tan u -|- tanh M cot w) 



(3.14) 



as already discussed, we will only focus on Eq. H3.13|l . 

To solve for 71, we will first impose a particular 0i, taking advantage of the fact that 0i can be thought of as a 
perturbation to the distant, Newtonian gravitational field of the black hole. We then integrate dji/dv with respect 
to v to find 7i(u,w), being careful to choose an appropriate integration constant by demanding that the perturbation 
vanish far from the black hole. The 71 (u, v) we find then automatically satisfies Eq. (|3.9() or (|3.14() . 
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In preparation for these calculations, it is useful to put this metric into Schwarzschild-like form. Plugging Eqs. 
and H3.12|l into Eq. H3.7|l . and making the coordinate transformation [29i| 

r = 2Af cosh^u/2, O^v (3.15) 

yields 

ds' = -e^^^ (^1 - ^) dp + e'^^-^^^ (^1 - ^) ' dr' + r'e'->^-'^^d0' + siiP Oe-^'^^dcjP . (3.16) 
Note that the transformation 13.15(1 imphes 



p = rsinOdl - — , (3.17) 

z = (r-M) cos 61 . (3.18) 

These relations can be particularly helpful when describing the perturbations in Schwarzschild coordinates. 

Although the line element l|3.16|l is technically exact, in the following analysis we will only work to first order in 
the perturbation. We thus should really expand e^"^^ ~ 1 + 2-0i, and likewise for 71. For compactness of notation, 
we leave these perturbations in the exponentials, with the caveat that they must be expanded to first order in all 
calculations. 

Before examining some interesting perturbations, it is useful to carefully study orbits in this general perturbed 
spacetime. Several results can be found in terms of ifii and 71, which facilitates later analysis. In the next section, 
we will examine the properties of equatorial orbits in this metric, and calculate the rate of periapse precession. 

IV. EQUATORIAL ORBITS OF BUMPY BLACK HOLES 

In this section, we study, as much as is possible, properties of bumpy black hole orbits that do not require specifying 
the perturbation. It is easy to find the geodesies in the equatorial plane {0 = 7r/2); for this first analysis, we will focus 
on this simple case. It is not, in fact, difficult to generalize the equations of motion to orbits beyond the equatorial 
plane. To solve these equations, however, appears challenging: the equations for r and 9 do not appear to separate 
in general (as they do, for example, in the Kerr case). We discuss this issue further in Sec. IVIIll 

As discussed in the Introduction, our goal is to understand how black hole bumpiness is imprinted upon measurable 
quantities — in particular, accumulated phases related to harmonics of the orbital frequencies. A simple (and 
historically important) effect is the shift of an orbit's periapsis. This shift is related to the mismatch between the 
r and (j) frequencies. In the weak-field limit, it describes perihelion precession, well-known from studies of planetary 
orbits in our own solar system. In this section, we will derive a differential formula for the periapsis shift. 

We begin by writing down a Lagrangian for orbiting bodies in this spacetime: put C = gabi'^i^ (where overdot 
denotes d/dr, with r proper time along an orbit). Note that C = —1. Since our focus here is equatorial orbits, we 
further put 6 — 7r/2 and 9 — 0. (Thanks to the reflection symmetry, we are guaranteed that an equatorial trajectory 
remains equatorial.) The result is 



e-^^'r^^ . (4.1) 



Varying C with respect to t yields 



d_ 

d^ 



e2V 



, (4.2) 



leading us to identify the orbital energy per unit rest mass ^ 



E = e^V-i ( 1 - ^ W (4.3) 



as a constant of motion. By varying with respect to 0, we likewise identify the orbital angular momentum orthogonal 
to the equatorial plane per unit ji: 



L = e-^'l''r^(t) . (4.4) 
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An equation for the radial motion follows from £ 



^2 ^ 



^2 _ g2^i 



2M 
r 



(4.5) 



It is easy to see that this reduces to the usual Schwarzschild equation of motion for V'l = 7i = [compare, e.g., Ref. 
|lfl| |. Eq. (25.16a)]. For the calculations we will perform momentarily, it is useful to multiply this by r^, defining 



R{t) 



^3.2 



(4.6) 



Using R{r), it is straightforward to remap the constants E and L to parameters that directly characterize the orbit. 
It is helpful to reparameterize r, 



P 



1 + e cos X 



(4.7) 



As X oscillates from to tt to 27r, the orbiting body moves from periapsis, rp = p/{l + e), to apoapsis, Va = p/{l — e), 
and back. In Newtonian gravity, this reparameterization facilitates identifying the orbits as closed ellipses; p is the 
orbit's semi-latus rectum and e its eccentricity. This intuition remains very useful in general relativity, though the 
ellipses do not close — the time for x to go from to 27r is greater than the time for the orbit to move through 27r 
radians of azimuth (j). 

The periapse and apoapse radii are the orbit's turning points. By definition, this means that f and hence R{r) 
equal zero at these points. This makes it possible to find E and L as functions of p and e: simultaneously solving 
R{rp) = and R{ra) = yields 



Eip,e) 



[(1 + e)2e2V'i^p) _ (1 _ e)2g2^iK)] [^2 _ + 4^/2(1 _ ^2)] 



p {e4V'(rp)(i + £:)2 _ 2M(1 + e)] - e4*('-")(l - e)^ [p - 2M(1 - e)]} ' 



(4.8) 



L{p, s) 



I e2^(rp) _ 2M(1 + £)] - e2'^('-») [p - 2Af (1 - e)] 

gii'ir^) (1 + £)2 _ 2M(1 + e)] - e4'A('-a) (1 - e)2 [p _ 2M(1 - e)] 



(4.9) 



For orbits that are circular (e = 0), rp = r^, so the conditions R{rp) — 0, R{ra) = do not provide separate 
information. In this case, we require dR/dr = 0; this condition guarantees that r remains zero. Solving R — Q, 
dR/dr — yields 



E{p, £ = 0) = e 



2M 
P 



p~2M~ p{p - 2M)dipi/dr 



p - 3M - 2p{p ~ 2M)di}}i/dr 



(4.10) 



L(p, £ = 0) = pe 



^^,ip) / M + p{p-2M)di;i/dr 

\ P-3M- 2p{p - 2M)d^i/dr 



(4.11) 



When the perturbation ipi is set to zero, these expressions correctly reduce to formulas characterizing E and L for 
Schwarzschild black holes: 



E{p,e) = 



/p2_4Mp + 4M2(l-£2) 
/ p[p-M(3 + £2)] 
1 - 2MIp 



v/1 - ^MIp 



(£ = 0) ; 



(4.12) 
(4.13) 



L{p,e) 



/ M 
p-M(3 + £2) 



/ M 
p~3M 



(£ = 0) 



(4.14) 
(4.15) 
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The function R{r) also determines the parameters of the last stable orbit: solving the equation 



dR 
dr 







(4.16) 



with E — E(j), e), L = L{p, e) determines a separatrix p(s) such that, at fixed eccentricity e, orbits with p < p(e) are 
dynamically unstable. (For circular orbits, we must solve d^R/dr^ — 0.) The result is rather messy, so we do not 
present it explicitly. When the perturbation is turned off, the separatrix is simply given by p = (6 + 2e)M . 

We are now ready to derive a general formula for bumpy black hole periapsis precession. First, note that, via Eq. 
1)4.7(1. the rate at which the orbital radius changes with respect to x is 



dr pe sin x 

dx (1 + ecosx)^ 

Combining Eqs. lll!5|l . and (|CT7|l . we find 



r s smx 
P 



(4.17) 



dcj) dcj) dr dr 
dx dr dr dx 

= e^'^i ['■(^^l L(p, e)e sin XP~ 



E{p,ef 



„2Vi[r(x)] 



2M 



.2^i[r(x)]L{p^ ef 



-1/2 



(4.18) 



This equation expresses, in differential form, the amount of azimuth that accumulates per unit x- Integrating d(j)/dx 
over < X ^ 27r gives the total accumulated azimuth in one full orbit. In Newtonian gravity, this number is 27r radians 
— Newtonian orbits are closed ellipses. The periapsis precession is the amount of "extra" azimuth that accumulates 
in one orbit: 



2ir 



dx^ - 2^ 
dx 



(4.19) 



Having derived about as many orbital properties as we can in general, much of the rest of this paper will be devoted to 
integrating this equation with particular choices of the perturbing functions ijji and 71 . We turn now to a calculation 
of these perturbations and an exploration of their effects. 



METRIC PERTURBATIONS I: POINT MASSES AT THE POLES 



In this and the following section, we calculate the first-order metric perturbation in two special cases: a pair of point 
masses at the poles (this section), and a ring of mass about the equator (following section). These two perturbations 
are particularly interesting because they produce a quadrupolar distortion of the spacetime. 

Point masses at the poles were analyzed by SPR in Ref. [l^. Their particular focus was on the perturbation at or 
near the hole's event horizon; we generalize their result to find the perturbation throughout the exterior spacetime. 

We begin by choosing, in Weyl coordinates (<, p, 2, 0), the ■0 perturbation 



m/2 



(5.1) 



As described previously, this choice has a simple interpretation: adding a point mass to the system changes the 
Newtonian potential at infinity — ip — by that of a point mass. Equation H5.1() describes a point with mass (as 
measured at infinity) /x/2 at z = 6, near the "north" pole. The corresponding tpi for a perturbing mass near the 
"south" pole is 



SP 



^P^ + {z + bf 



(5.2) 



The complete tp perturbation is given by adding the contributions from the north and south poles. 

Our choice forj^ has a curious property: it corresponds to a nonspherical naked singularity at p — 0, z = ±b (the 
Curzon solution |33j|). One might object that a spacetime which includes such a singularity can not be physical, since 
naked singularities are presumed, by cosmic censorship, not to exist in nature. Therefore, perhaps we should exclude 
such spacetimes from consideration. 
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As we have repeatedly emphasized, given the goals of this analysis, the physicality of our trial spacetime is irrelevant. 
We seek a family of spacetimes that we can compare with those found in nature. In particular, we ned a family of 
spacetimes that deviate controllably and only slightly from black holes. For the purposes of uncovering whether a 
spacetime is a black hole or not, it makes no difference whether the deviations come from physics we expect or not. 
We merely want to set limits on the extent of any deviations that might exist. 

Having specified ipi, we turn to the computation of 71. Focus first on the perturbation at the north pole. Switching 
to the prolate, spheroidal coordinates u and v, we rewrite 



v^r = - 



m/2 



[iVP sinh^ u sin^ v + {b — M cosh u cos v)^] 



1/2 



(5.3) 



Insert this into Eq. (|3.13|l : 



97? 



NP 



dv 



/iM (6 — M cosh u cos v) sin v 
[AP sinh^ u sin^ v + {b — M cosh m cos i;)^] '^^'^ 
fi (6 — coshucosw) sinu 

M r . , 2 ~2 7t T „l3/2 " 

smh Msm u + (0 — cosh w cos wj^ 



(5.4) 



We have defined b = b/AI. 

We now integrate this with respect to v to compute 7^^^. Several manipulations help. First, simplify the denomi- 
nator: 

sinh^ wsin'^ w + (6 — coshucosw)^ = sinh^ m + 6^ — 26 cosh u cos w + cos^ w 

= ( cos V — b cosh uj +b — b cosh u + sinh u 

= ^cos u — 6 cosh + ^1 — sinh^ M . (5-5) 

Next, put X = cos V, so dx ~ — sin v dv: 

I —6 + a; cosh u] dx 



X — b cosh u 



„,NP M 

^1 =M 



Let y — X — 6 cosh u; dy = dx (for these purposes u is fixed): 



1 — 6^ 1 sinh^ u 



3/2 



(5.6) 



„,NP 



JJ_ 
M 



(hsi 



sinh u + y cosh u J dy 



?/2 + ( 1 — 52 j sinh^ u 
y cosh u dy 



3/2 



5 sinh u dy 



-I- ^1 — b'^^ sinh^ u y'^ + (^1 — 5^^ sinh'^ u 



3/2 



(5.7) 



The first of these integrals is: 



y cosh u dy 



cosh? 



y2 + ( 1 — 52 j sinh^ u 



3/2 



1/2 + (1 — fe2^ sinh^ u 
coshu 



1/2 



sinh u sin i; + ( 6 — cosh u cos 1 



1/2 • 



(5.8) 
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The second yields 



b sinh u dy 



yb sinh u 



1-52 sinh 



3/2 



1 — 52 I sinh u 
b 



1 — 6^ I sinh 



1/2 



cos V — b cosh ? 



1-52 



sinh^ u sin^ 



6 — cosh u cos w 



1/2 



(5.9) 



Restoring factors of M, the complete integral is 



fiM {b cos V — M cosh u) 



(Af 2 — 52 -J j^,j2 sjjj}^2 ^ gj^2 „ ^ _ cosh u cos ti)^ 



1/2 



(5.10) 



In these integrals, we have neglected integration constants that must now be determined. We do so by requiring 
that 7^^^ go to zero at large radius, i.e. as w — s- cx). Our current form of the solution has 



iJLM 



hp - 62 



(5.11) 



This is the value that must be subtracted to guarantee that 7^^ is well behaved at large radius. We thus finally 
obtain 



M^-b^ 



b cos V — M cosh 1 



M2 sinh^ u sin^ v + (b — M cosh u cos vY 



1/2 Af2_52 



(5.12) 



To account for the south pole perturbation, we just add the perturbation with b —b. [Note that we do not add the 
south perturbation with an overall minus sign, as is done in Ref. [2^. That sign choice only holds for u ^ 1; see Ref. 
^9], discussion following Eq. (3.14). Note that this has no impact on our ability to study strong field orbits — the 
region u ^ 1 corresponds to a region just barely outside the event horizon, well inside of the last stable orbit.] By 
direct substitution, one can verify that this solution for 71 satisfies Eq. (13 .911 as well. 

The perturbed metric is now fully specified. For our purposes, it is very convenient to use Eq. H3.15|l to convert to 
Schwarzschild coordinates. Doing so and defining the function 



D(r, 6*, b) = (r2 - 2Mr + b^ + IVP cos^ 9 - 2br cos 9 + 2bM cos 9) 
the complete perturbations V'l and 71 are 



1/2 



^1 



2D{r,9,b) 2D{r,9,^b) 



(5.13) 
(5.14) 



71 - 



A-p - 62 



M + 6cos6'-r\ f M -bcos9 ■ 



Dir,9,b) J V D{r,9,~b) 



2iiM 

^ (5.15) 



M2 - 62 



These expressions hold for all r > 2M . Note the Schwarzschild coordinate locations of the point masses are r = b + M, 
9 — 0,Tr. We must choose b > M in order for the perturbing masses to be in the hole's exterior. 

Combining Eqs. (|4.8|l . H4.9|) . (|4.18|) . H4.19|l . H5.14|l . and H5.15|l . we can now calculate the periapse precession for 
equatorial orbits and the two-point-mass perturbation. Numerical results for the general, strong-field case will be 
discussed in Sec. IVII Bl here, we show the weak- field result {p S> M, p b). This weak-field phase shift breaks 
naturally into 3 pieces: 

^^points ^ A</.Schw(M + /i) + A(/)anom(M) + A'/'prollM, &) , (5-16) 

where 

A0schw(A/ + /i) = \ —5 (18 -fe) , (5.17) 

p Ap'^ 

A0anom(A^) = -^(l + 2£2) , (5.18) 



A0proi(/x,6) = (5.19) 
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These three pieces each have a simple, physical explanation. The first, A^schw is just the periapsis precession 
expected for a particle orbiting a Schwarzschild black hole with total mass M + fj, (to leading order in /z, and to 
second order in 1/p). The second, Ai^anom, is an "anomalous" contribution to the precession arising from the fact 
that our perturbation is non-spherical and has non-zero total mass. As we show later, this term can be eliminated by 
appropriately designing the perturbing mass distribution. 

The third piece, A(/)proi, is the leading order contribution to the precession that arises from a prolate, quadrupolar 
distortion to the black hole. This term agrees exactly with a Newtonian calculation of periapse precession in the 
presence of a prolate quadrupolar distortion (cf. Appendix IXjl. 



VI. METRIC PERTURBATIONS II: AN EQUATORIAL RING 



We next examine perturbations due to a circular, equatorial ring of mass p about a black hole. Our procedure is 
essentially the same as that used in the two-point-mass example, though the form of the perturbing potential makes 
the calculation a bit more complicated. We choose for the potential of a ring with radius p = b and total mass at 
infinity p.: 



2Tr 



+ 2;2 + 52 „ 2bp COS ^] 



1/2 



(6.1) 



This potential can be re-expressed using the complete elliptic integral of the 1st kind; it turns out to be more convenient 
to leave it as written here. 

Inserting this into Eq. H3.13|l and integrating yields 



7i 



g{u,v,£,) h{u,v,^) 



where 



/(u,w,^) — /icoshu (6^ -I- il/^ cosh'^ It — M6sinhusinucos^) , 
g(u,w,0 = TT [(M^-l-b^cos^Cjcosh^u-f fe^sin^^] , 



[M sinh usin v + M cosh mcos v + b — 2A/6 sinh u sin u cos 



1/2 



(6.2) 

(6.3) 
(6.4) 

(6.5) 



The integration "constant" C(^) is chosen, as in the two-point-mass case, to make 71 vanish at large distances. It has 
the value 



pM 



TT (Af2 + 62 C0S2 



(6.6) 



To evaluate the periapsis precession, it is useful to have these results in Schwarzschild coordinates. Using Eq. 13.15|1 . 
/, 5, and h become 



f{r,6,^) = — (r - M) \b^ + (r - M)^ - b\/r^ - 2Mrsm9cos^ 



1 + cos^ e ) (r - + &2 sin^ ^ 



(r^ - 2 A/r) sin^ + {r - Mf cos^ 9 + b^ ~ 2bVr^ - 2Mr sin 6 cos ^ 



1/2 



(6.7) 
(6.8) 
(6.9) 



The ring's radius in Schwarzschild coordinates is r = M + \J Af 2 + 52 ; any choice 6 > will produce a ring in the 
hole's exterior. 

Since the integration over ^ commutes with other operations (notably expanding for large p) it is straightforward 
to combine Eqs. (|4.8() . (|4.9() . I|4.18|l . H4.19() . (|6.1|) . and H6.2|l to compute the periapsis precession in this spacetime. We 
do so numerically in Sec. IVllBl for now, we evaluate the weak- field precession [p ^ Af, b): 



A0""s = A0Schw(M + A') + A0 anom 

(p) + A(t)oh\{p,b) 



(6.10) 



The terms A0schw(A^ + Ai) and A(f)anom{p) are defined in Eqs. ()5.17fl and ()5.18|l : they are identical to the corresponding 
terms for the point mass perturbations. The term 



A(j)ohiip,b) 



' 2Afp2 



(6.11) 
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gives the contribution to the precession due to an oblate quadrupolar distortion to the black hole. As in the prolate 
case, this term agrees exactly with a Newtonian calculation of periapse precession presented in Appendix fXl 

VII. METRIC PERTURBATIONS III: PURE QUADRUPOLE PERTURBATION 

A. Weak-field analysis 

Each of our perturbations changed not only the quadrupole {I = 2) structure of the spacetime, but the spacetime's 
total mass {I — 0) as well. To construct a perturbation that is purely quadrupolar (at least in the weak field), we take 
advantage of the linearity of the perturbed Einstein equations to superpose a ring of mass ±/i with a pair of points 
of mass TA'/S each. A negative mass perturbation may seem strikingly unphysical, but as discussed in Sec. IIIII our 
goal is not to build perturbations that are likely to exist in nature. (Given the naked singularity interpretation, it is 
arguable that a negative mass perturbation is no less physical than one of positive mass.) 

From now on, we will focus on this zero-mass perturbation. The weak field result periapse precession for this case 

is 

A0b°th ^ A0g^,^„(M) + A0quad(M, b) , (7.1) 

where the quadrupole contribution is 

A<^,,.d(A^,5)=±g^. (7.2) 

Notice that the "anomalous" contribution to the precession does not appear in this case. Making the perturbation 
have zero mass apparently suffices to eliminate that term. What remains is a periapsis precession term corresponding 
to a purely quadrupolar deformation of the spacetime. We have succeeded in constructing a spacetime that is almost 
that of a black hole, that is good deep into the strong field, but has one multipole moment — M2 — with the "wrong" 
value. 

B. Strong- field analysis 

Having demonstrated that our calculations match up in a sensible way with Newtonian and relativistic calculations 
of periapse precession, we now show numerical calculations for the strong fields of bumpy black holes. In all of our 
figures and discussion, our focus is on the precession induced relative to periapsis precession for a normal black hole. 
To help calibrate the effect that the bumpiness has, we show in Fig. ^ the periapsis precession (50 (the amount of 
azimuth (j) accumulated in one orbit, minus 2tt) in the strong field of a Schwarzschild black hole. Notice that S(j) 
becomes quite large in the strong field, corresponding to almost 2 extra revolutions around the hole. 

Figures [3 and 13 show the shift Acj) associated with two different choices of bumpiness. In both plots, we have put 
fi — 0.01; since the Einstein equations and the equations of motion all scale linearly with /i, A4> is likewise linear in 
fi. Thus, it is a simple matter to rescale to other masses. We examine strong-field precession for e = 0.7 and s = 0.1 
in these plots; we also compare these results to the weak-field prediction (|7.2|) . 

Several features are evident. First, in both cases, the strong-field results asymptotically approach the weak-field 
formula for large p — an important sanity check. Interestingly, as b is increased from 1.5A/ to 2.5Af, the weak- 
field formula changes from overestimating the periapsis shift to underestimating it. This behavior may be due to 
deviations of our potential from that of a perfect quadrupole. Second, at most orbits, the eccentricity has very 
little impact on this shift, at least for the values we examine. This is not surprising; the bare Schwarzschild periapsis 
precession (Fig.^l is likewise fairly insensitive to the eccentricity. Finally, even into the strong field, the periapsis shift 
scales approximately with b^ (the weak-field prediction). We can thus regard ^b'^/M^ as a measurable, dimensionless 
"bumpiness parameter" for the black hole. 

Most interesting is the robustness with which the quadrupolar bump is manifested in the periapsis result: the effect 
is quite pronounced in the strong field. Indeed, the weak-field prediction (|7.2|) underestimates the degree of precession 
due to the hole's bumpiness by a factor ^ 3 over much of the strong field (~ 10 as we approach the last stable orbits). 
Over much of the strong field, the effect is large enough that, by observing over multiple orbital cycles, we should be 
able to set very interesting limits on the bumpiness of black hole candidates. For example, at p ~ 20Af , the periapsis 
shift per orbit is roughly 
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A bumpiness ^b'^/M^ ~ 0.02 should thus be easily measurable after tracking roughly 1000 orbits (corresponding to 
when the bump shifts the accumulated phase by about 1 radian). Conversely, tracking the phase for N orbital cycles 
at this value of p should make it possible to constrain the bumpiness to be 



Better limits can be obtained for orbits at smaller p. 

Obviously, more detailed work is needed to carefully examine how well black hole bumpiness can be measured in a 
variety of scenarios (akin to Ryan's analysis of multipole measurability by LISA, Ref. 26]). But, these results suggest 
that measurements which coherently follow orbital phases — such as gravitational-wave measurements and x-ray 
timing — should be able to place stringent constraints on the bumpiness of black hole candidates. Bumpy black holes 
should thus be very useful tools in designing a formalism to map the strong field structure of black hole candidates 
in nature. 



In this paper, we have laid the foundations for a null experiment to test whether a massive compact object is a 
black hole. The bumpy black hole spacetimes we construct differ only slightly from normal black hole spacetimes; 
and, the difference is controlled by a simple adjustable parameter — the hole's "bumpiness" . It should be possible to 
compare the properties of black hole candidates in nature with these bumpy black hole spacetimes. If these objects 
are in fact the black holes of general relativity, measurements will show that the natural spacetimes have a bumpiness 
of zero. 

Quite a bit more work is needed in order to make the bumpy black hole concept useful in practice for astrophysical 
measurements: 

• Foremost is the need to generalize this analysis to bumpy Kerr black holes — zero angular momentum is a highly 
unrealistic idealization. In Appendix 151 we show how, by choosing tp, 7, and using an appropriate coordinate 
transformation, the Weyl metric H3.1|l encompasses Kerr black holes (in Boyer-Lindquist coordinates). There 
should then be no severe difficulty perturbing this metric to build bumpy Kerr black holes, though the details 
are likely to be complicated. 

• Probably next in importance is generalizing the orbits which we analyze to include inclination with respect to 
the hole's equatorial plane. Besides the astrophysical motivation — we do not often expect orbits to be confined 
to a special plane — the inclusion of an extra degree of orbital freedom offers opportunity. Motions out of the 
plane are characterized by oscillations with a frequency ilg which is generically different from the frequencies 
and fir discussed in this paper. These oscillations thus offer an additional "handle" by which deviations from 
the black hole spacetime can be characterized. 

As already mentioned in Sec. II VI the equations of motion for inclined orbits of bumpy black holes do not appear 
to separate (as they do for normal Kerr black holes). However, the fact that, by definition, bumpy black hole 
spacetimes are nearly the spacetimes of black holes suggests that the equations of motion must nearly separate. 
In other words, the degree to which the 9 motion couples to the r and 4> motion must be small — no doubt 
controlled by a coupling factor that is of order n/M. It may be possible to take advantage of this smallncss to 
usefully describe inclined and eccentric bumpy black hole orbits. 

• We have focused on perturbing mass distributions that produce, in the weak field, a purely quadrupolar spacetime 
distortion. We chose to focus on this case because an incorrect quadrupole moment is sufficient to falsify the 
black hole hypothesis. There is no reason why we could not go beyond this: by using rings out of the hole's 
equatorial plane, one could imagine building essentially any multipolar distribution whatsoever. Indeed, the 
fact that the equation for the metric function -0 in Weyl coordinates [Eq. (|3.2|l ] is simply V^i/; = tells us that 
it is simple in principle to specify perturbations whose weak-field multipolar structure is completely arbitrary: 
the perturbation 



will work perfectly. The parameter Bi is a generalized Z-polar "bumpiness" . With this ansatz to define our 
deviations, we can build a bumpy black holes with almost arbitrarily shaped bumps. This will make it possible 
to strongly constrain the properties of black holes in nature. 



VIII. CONCLUSION 




(8.1) 
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As this paper was being completed, an analysis appeared on the gr-qc archive by Ashtekar et al. [35| of the 
multipole moments of isolated horizons |36j . Although we have not investigated this in any depth, it may be 
beneficial to pursue a connection between the bumpiness of a black hole and the multipole moments expressed 
in the language of Ref . [s^ ■ 

With these generalizations in hand, it should be possible to begin examining detailed mechanisms by which orbital 
frequencies can be imprinted on astrophysical observables. For example, one can imagine analyzing accretion disk 
models to see how a spacetime's bumpiness is imprinted on quasi-periodic oscillations in a source's x-ray spectrum 
[T^ IstI Is^ . Another example is in gravitational- wave science. For these ideas to be useful for testing the nature of 
black hole candidates, we will need to model the gravitational-wave emission and inspiral of compact bodies captured 
by bumpy black holes. This problem may not be much more difficult than the corresponding problem for normal 
Kerr black holes — because the wave amplitude is itself perturbatively small, the inspiral and wave generation should 
decouple (at least to first order) from the spacetime's bumps. We hope to address at least some of these issues in 
future work. 
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APPENDIX A: PERIAPSE PRECESSION IN NEWTONIAN THEORY 

An elliptical orbit precesses even in Newtonian gravity if the orbit is about a body with a quadrupole moment. 
Here, we calculate this precession, presenting results in a form useful for making contact with this paper's relativistic 
results. 

A body with a mass and a quadrupole moment has a Newtonian gravitational potential 

^_M_32^^ (Al) 

For further discussion, see Thome's voluminous treatise on multipole moments in general relativity |34|. The tensor 
lab is the symmetric, trace-free quadrupole moment of the gravitating source: 

1ab = J d^r p (^XaXt- ^Sabr'^^ ■ (A2) 

The vector n"" is a direction cosine. For notational simplicity, we put lab'n'^n^ = Q- 
A body in an equatorial orbit around this object obeys a simple equation of motion: 

drV 

— = E-^ 

dt / 2r2 



Rncwtir) 



(A3) 



where L is the component of orbital angular momentum perpendicular to the equatorial plane. We are interested in 
eccentric orbits, so we reparameterize in the usual way: 

P (A4) 



1 + e cos X 
which is equivalent to 

w = i;(l H- ecosx) (A5) 

with u = 1/r, V ^ l/p. The turning points of an eccentric orbit are at apoapsis, Ua — v{l ~ e), and periapsis, 
Up = v{l + e). Solving the equation of motion at both turning points — Rnewt{ua) = 0, i?ncwt(wp) = — yields a 
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solution for the energy E and angular momentum L as functions of v and e: 

E = l{e^-l)Mv + l{e'-l)'Qv\ (A6) 
L2 = ^+3^3 + ,2)Q„. (A7) 

V Z 

An orbit corresponds to moving through the range x = to x = 27r. We want to find the amount of azimuthal angle 
(j) that accumulates over that orbit; that number (minus 27r) is the accumulated periapsis precession. First, use the 
above results for E and L in the radial equation of motion: 

= ie2„sin2x[2M + 3(£^-3)Qt;2-6eQf2cosx] (A8) 
Next, using r = l/v{l + £COSx) we have 

dr ^ . dx , , „N 

_^,,,smx^. (A9) 



We obtain 

'dxV 2M + 3{e^ ~3)Qv -GeQvcosx 



dt J 2vr^ 



(AlO) 



We leave the in the denominator to cancel another factor that will appear shortly. 

To connect this to the azimuthal angle (/), we use the fact that L = r"^ dip/dt. Combining this with (|A10p . we obtain 
an equation for the differential periapsis advance: 



^ ^ d(f>\^ f dx^ ^ 



dxj \dt J \dt 



L 



2M + 3(e2 _ 3)Qw - QeQvcosx 
2M + 3(3 + £2)gt,2 



2M + 3(e2 - 3)Qi;2 - GeQv"^ cosx 



(All) 



We now take the square root. We assume that Qv^ = Q/p^ ^ M and expand in v. This amounts to a weak field 
expansion. We find, to leading order in Qv^ /M , 



d0 ^ 3Qt;2(3 + £cosx) 

dx 2M ^ ' 



Integrating over an orbit yields 



2ir 



We now compute Q for the cases examined in the text: two point masses on the symmetry axis, and an equatorial 
ring. The direction cosines are easy to calculate; we are only interested in orbits on the equatorial plane, so = 0. 
By axisymmetry, all directions in the equatorial plane are equivalent, so we choose = 1, = 0. Then, Q = Txx- 

For a point mass /i/2 at the north pole, 

ir = ^diag(-l,-l,2); (AM) 

for a point mass at the south pole, we find the same result. The pair of point masses thus has Q = Using 
this Q in Eq. I|A13|I we find 

^^points ^ _37rp^ _ 
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Repeating this exercise for the ring, we find 



so Q = ^J■b^/6, and 



i;r = ^<ii«s(i.i,-2), (A16) 



A*'™ = - (A17) 



Notice that A(/)P°'"''' = — 2A0'''"s. This follows from the spherical harmonics that describe the ring [F2o(7''/2)] and 
the point masses [l2o(0) = -~2Y2o{tt/2)]. 

APPENDIX B: THE KERR METRIC FROM THE WEYL METRIC 

In this appendix, we lay out the coordinate transformations and the choices of ip and 7 needed to go from the Weyl 
metric to the Kerr metric in Boyer-Lindquist coordinates. Although only the uncharged version is astrophysically 
relevant, we show the results for general charge Q. This calculation shares much with the Newman- Janis algorithm 

We begin with the Weyl metric in prolate spheroidal coordinates, Eq. H3.7|) . We then choose 



Af coshu = r - M + i^/a^ + Q2 . (Bl) 

Note that u must be complex; the black hole parameters M, a, and Q, as well as the radial coordinate r are of course 
real. As we discuss further below, this implies a severe constraint on the values u may take in the complex plane. 
It follows that 

dr ^ M sinh u du , 
M^cosh^u EE M^|coshu|^ 

- (r-M)2 + a2 + Q2^ 
sinh^ u = - 2Mr + + . (B2) 

We also choose v = 9. For notational convenience, define 

A = - 2Mr + + Q'^ , 

E = + cos^ 6 . (B3) 

Substituting into the metric H3.7() . we find 



(B4) 



Next, choose 

A 



g2V ^ 



e'"" = ^ T- , (B5) 

A + M2 sin^ 61 ' ^ ^ 

and put 

= dt' - asm^0d(t)', 
d(f, = (ll±^ dcP' -^dt' . (B6) 

Dropping the primes on t and 4>, we see that the Weyl metric reduces to the Kerr-Newman metric in Boyer-Lindquist 
coordinates: 

ds^ = {dt - a sin^ 9 dcp)^ + [(r^ + a^) dc/) - a dtf + ^ dr^ + dO^ . (B7) 

Zj Za 



19 



Compare, for example Ref. [Ta|, Eq. (33.2) and Ref. W\, Eq. (12.3.1). 
In prolate spheroidal coordinates, the functions ijj and 7 are given by 



r-^(u) + cos^ v 
M2 cosh^ u + M2 sin^ u ' 



where 



r{u) = M (cosh M - 1) - i + Q2 _ (BIO) 

The requirement that r(u) be real implies a severe constraint between the real and imaginary parts of u: we must 
have Im(A/coshu) ~ +i\/a'^ + Q^. Writing u ~ + iui, we find 



Msin (uj) sinh (ur) = \/ a'^ + . (Bll) 

In other words, u cannot take just any value in the complex plane, but must be confined to the trajectory defined by 
Eq. ||BTT|| . 

Using these results, it should straightforward to generalize our calculations to describe bumpy Kerr black holes. 
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FIG. 2: Periapsis precession in the strong field of a bumpy black hole: the "extra" azimuth accumulated due to the hole's 
bumpincss as a function of {p,e). Wo have put fi — 0.01 and b = 2.5. By the linearity of all relevant equations in mass, the 
degree of precession scales proportional to jj,. The solid curve is for e = 0.7; the dashed one is for e = 0.1. The dotted curve 
shows the weak-field prediction. Both strong-field results asymptotically approach the weak-field prediction for large p. 
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FIG. 3; Periapsis precession in the strong field of a bumpy black hole: the "extra" azimuth accumulated due to the hole's 
bumpiness as a function of (p,e). This plot is identical to Fig. |21except that we have put b = 1.5M. 



